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1. Introduction 

The proton- proton fusion process, pp — > de+fg, is a fundamental reaction for the nuclear 
astrophysics, especially important for the understanding of the star evolutions [1] and solar 
neutrinos [2, 3, 4]. However, the process has never been studied experimentally because 
the event is extremely unlikely to take place in the laboratory at the proton energies in 
the sun. The calculation of the transition rate and its uncertainty has naturally become a 
challenge to nuclear theory. The first calculation of the process was carried out by Bethe 
and Critchfield [5] in 1938. This estimation was improved by Salpeter [6]^ in 1952. Later, 
small corrections, such as the electromagnetic radiative corrections, were considered by 
Bahcall and his collaborators [8, 9] in the framework of effective range theory. Recently, 
accurate phenomenological potential models were employed to study the process [10, 11]. 
Furthermore, in Ref. [12] the two-nucleon current operators were calculated from heavy- 
baryon chiral perturbation theory (HBxPT) up to next-to-next-to-next-to leading order 
(N'^LO), and Park et al. obtained quite an accurate estimation (~ 0.3% uncertaity) 
for the process by fixing an unknown parameter, so-called low energy constant (LEG), 
which appears in the two-nucleon-axial-current contact interaction in terms of the tritium 
lifetime [13, 14]. 

The kinetic energy relevant to the pp fusion process at the core of the sun is quite low, 
kTc — 1.18 keV, where Tc is the core temperature of the sun, Tc — 13.7 x 10^ K, and 
k is the Boltzmann constant. The proton momentum at the core, Pc — ^2mpkTc ~ 1.5 
MeV, where Trip is the proton mass, is still significantly small compared to the pion mass, 
m^r — 140 MeV. Therefore, we may regard the pion as a heavy degree of freedom for the 
pp fusion process. It may be convenient and suitable to employ a pionless effective field 
theory (EFT) [15], in which the pions are integrated out of the effective Lagrangian for the 
process in question. The pp fusion process in the pionless theory has been studied by Kong 
and Ravndal [16] up to next-to leading order (NLO) and by Butler and Chen [17] up to 
fifth order (N^LO). Thanks to the perturbative scheme in EFT, the accuracy of the N^LO 
calculation would, in principle, be (Q/A)^ ~ (1/3)^ — 1%, where Q/A ~ 1/3 is a typical 
expansion parameter in the pionless theory. However, because of lack of the experimental 
data to fix an unknown LEG Lia which appears in the two-nucleon-axial-current contact 
interaction in the pionless effective Lagrangian, an uncertainty estimated in the pionless 
EFT for the pp fusion process is still significantly larger than what is expected from the 
counting rules of the theory. 

In this work, we employ a pionless EFT with di-baryon fields [18, 19, 20]. The 
amplitude for the pp fusion process at the zero proton momentum is calculated up to 
NLO. We introduce two di-baryon fields [24], which have the same quantum numbers 
as those of S-wme two-nucleon states (^5*0 and ^5*1 states), as auxiliary fields: after 
integrating out the di-baryon fields we do have the ordinary pionless theory without the 
di-baryon fields. However, as have intensively been discussed in Refs. [18, 19, 25, 26], 

^For a recent historical review, see Ref. [7] . 

■^Wc have employed the same formalism in the studies of the two-body processes, such as neutron- 
neutron fusion [21], radiative neutron capture on a proton at BBN energies [22], and neutral pion pro- 
duction in proton- proton collision near threshold [23] . 



2 



with the aid of the di-baryon fields, resummation of the effective range correction up 
to the infinite order is naturally introduced, which greatly simplifies the calculation of 
higher order corrections to the wave functions. In addition, the new counting rules make 
the expansion parameter Q much improved, and it is not necessary to employ the power 
divergence subtraction scheme [27] any longer. Furthermore, by assuming that the leading 
order (LO) contribution in the di-baryon-di-baryon-current contact interaction can be 
determined mainly from the one-body current interaction as discussed in Ref. [19], we can 
reproduce the results from the effective range theory [28] in the LO calculations of the 
pionless EFT with the di-baryon fields. The NLO correction, the di-baryon-di-baryon- 
current contact interaction denoted by the unknown LEG Iia, is approximately presumed 
to be the two-body (2B) current correction in the pionful calculations. We fix the LEG 
liA by using the relative strength of the two-body matrix element to that of the one-body 
contribution, S2B [14], which has been determined from the accurate tritium lifetime 
datum. (We discuss it in detail later.) Gonsequently we can make our estimation of the 
pp fusion amplitude free from unknown parameters. Moreover, though our calculation is 
rather simple and is only up to NLO, we can obtain a result comparable to that from the 
accurate potential model calculation within the accuracy better than ~ 1%. 

This paper is organized in the foUowings: in Sec. 2, we introduce the pionless effective 
Lagrangian with the di-baryon fields up to NLO, and in Sec. 3, we fix the LEGs which 
appear in the initial and final two-nucleon states by using the effective range parameters. 
In Sec. 4, the amplitude for the pp fusion process is calculated up to NLO. We show our 
numerical results in Sec. 5. In Sec. 6, discussion and conclusions are given. 

2. Pionless effective Lagrangian with di-baryon fields 

For the low-energy process, the weak-interaction Hamiltonian can be taken to be 

n = ^i.J", (1) 

where Gp is the Fermi constant and Vud is the CKM matrix element, is the lepton 
current 1^ — UeJij,{l — "y5)vu, and is the hadronic current. We will calculate the two- 
body hadronic current J'^ from the pionless effective Lagrangian with di-baryon fields up 
to NLO. 

We adopt the standard counting rules of pionless EFT with di-baryon fields [18]. In- 
troducing an expansion scale Q < A(~ rriT^), we count the magnitude of spatial part of 
the external and loop momenta, \p\ and \l\, as Q, and their time components, p° and 
as Q^. The nucleon and di-baryon propagators are of Q^^, and a loop integral carries 
Q^. The scattering lengths and effective ranges are counted as Q ~ {7, 1/ao, 1/pd, ^/f^o} 
where 7, ao, Pd and vq are the effective range parameters for the 5'-wave A^A^ scattering; 
7 = y/rrij^B, where B is the deuteron binding energy, qq is the scattering length in the ^Sq 
channel, pd and Tq are the effective ranges in the ^Si and ^5*0 channel, respectively. The 
orders of vertices and transition amplitudes are easily obtained by counting the numbers 
of these factors in the Lagrangian and diagrams, respectively. As discussed below, some 
vertices acquire factors like Tq and pd after renormalization and thus their orders can differ 
from what the above naive dimensional analysis suggests. Note that we do not include 
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the higher order radiative corrections, such as the vacuum polarization effect [29] and the 
radiative corrections from one-body part [30]. 

A pionless effective Lagrangian with di-baryon fields may be written as [18, 19] 



St 1 



(2) 



where Cn is a one-nucleon Lagrangian, Cg is the spin-singlet {^Sq state) di-baryon La- 
grangian including coupling to the two-nucleon, Ct is the spin-triplet {^Si state) di-baryon 
Lagrangian including coupling to the two-nucleon and Cgt describes the weak-interaction 
transition (due to the axial current) from the ^Sq di-baryon to the di-baryon. 
A pionless one-nucleon Lagrangian in the heavy-baryon formalism reads 



= NUiv ■ D -2igAS ■ A + 



1 



2m 



N 



[v ■ Df 



+ 



N, 



(3) 



where the ellipsis represents terms that do not appear in this calculation, v'^ is the velocity 
vector satisfying v'^ = 1; we choose f'* = (1, 0), and 5"^ is the spin operator 25''* = (0, a). 



Covariant derivative reads as i?^ 



d - -T 

2 ' 



where is the external isovector 



vector current, and = ~2^' where A/j, is the external isovector axial current, qa 
is the axial- vector coupling constant and mjv is the nucleon mass. 
The Lagrangians that involve the di-baryon fields are given by 



^st — 
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( ny+pd\ 



1 
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4(Ar^Pi'5°)Ar) + h.c.l , (4) 
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U-yt\t\{N'^Pf^'^N) + \,.c. 
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lA 



sluA^ + h.c. 



(5) 
(6) 



where Sa and ti are the di-baryon fields for the ^5*0 and ^5*1 channel, respectively. The 
covariant derivative for the di-baryon field is given by = dn — iCV^^^ where V^^^ is 
the external vector field. C is the charge operator for the di-baryon field; C — 0, 1, 2 
for the nn, np, pp channel, respectively. as,t is the sign factor ag^t = ±1 and A^ ^ is the 
mass difference between the di-baryon and two nucleons, m^^t = 2mAr -l- A^^^. ys^t is the 



di-baryon-two-nucleon coupling constant. Pj 
or channel; 



is the projection operator for the S — ^Sq 



1 

78 



p.' 



1 

V8 



(7) 



where (Tj (to) is the spin (isospin) operator. Note that, as mentioned in the Introduction, 
we separate the di-baryon-di-baryon-current contact interaction in Eq. (6) into the LO 
and NLO terms. The LO interaction proportional to qa is determined by the one-body 
axial-current interaction and the factor ^{tq + Pd)/ ^/roPd is included so as to reproduce 
the result from the effective range theory at LO. The NLO correction is parameterized 



4 




Figure 1: Diagrams for the dressed di-baryon propagator including the Coulomb inter- 
action. A double-line with a filled circle denotes the renormalized dressed di-baryon 
propagator. Double-lines without the filled circle and single-curves denote the bare di- 
baryon propagators and micleon propagators, respectively. Two-nucleon propagator with 
a shaded blob denotes the Green's function including the Coulomb potential. A (spin- 
singlet) di-baryon- nucleon-nucleon (sNN) vertex is proportional to the LEC y^. 



by the LEC Iia- More detailed discussion about the separation of LO and NLO contact 
interaction with external probe in the di-baryon formalism can be found in Ref. [19]. 

3. Initial and final A^A^ channels 

The typical energy of the pp fusion reaction is very low, as discussed in the Introduction, 
so we can assume that the dominant channel of the reaction is from the initial ^5*0 pp 
state to the final ^Si deuteron state. In this section, we fix the LECs which appear in 
the initial and final two-nucleon states for the pp fusion process from the effective range 
parameters. 

In Fig. 1, LO diagrams for the initial pp state in ^5*0 channel, i.e., the dressed ^5*0 chan- 
nel di-baryon propagator, are shown where the two-nucleon bubble diagrams including the 
Coulomb interaction are summed up to the infinite order. The inverse of the propagator 
in the center of mass (CM) frame is thus obtained by 

iD;\p) = ias{E + 5s)-iylJo{p), (8) 

with 

where is the outgoing two-nucleon Green's function including the Coulomb potential, 

G^c\e)- ^ f,\y ^. ■ (10) 
E - HQ-Vc + ie 

E is the total CM energy, E ~ p^/rriM. Hq is the free Hamiltonian for two-proton. 
Ho = p'^/niN, and Vc is the repulsive Coulomb force Vc = a/r: a is the fine structure 
constant. Employing the dimensional regularization in d = 4 — 2e space-time dimension, 
we obtain [31, 32] 



Mp) = 



Stt 



N , 



-^Kv)-c,—{ip), (11) 
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Figure 2: Diagram for the S-wscve pp scattering amplitude with the Coulomb and strong 
interactions. See the caption of Fig. 1 for details. 



where jj, is the scale of the dimensional regularization, Ce = 0.5772 • • •, and 

1 



h{rj) — Re il>{irj) — larj , Re il>{rj) — rj^^2 



— Ce , 



2 ^ 27ir] 

y g27rr) _ I ' 



arriN 



Thus the inverse of renormalized dressed di-baryon propagator is obtained as 



At: 



■p + aniNiiyr]) + ipC^ 



mNVs 



1 9^ 



where is the renormalized mass difference 



am 
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Stt 



- - SC^; + 2 + In 
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TTfl 



(12) 
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In Fig. 2, a diagram of the 5- wave pp scattering amplitude with the Coulomb and 
strong interactions is shown. Thus we have the S-wave scattering amplitude as 



{-Ws^o)iDs{p){-iys^o) 
. Aw 



TUN _ _ amNh{ri) -ipC^' 



(15) 



mNVs 



with 



/ 



(2^ 



(+)\ 



(2^ 



(16) 



where {k\il)^^'') are the Coulomb wave functions obtained by solving the Schrodinger equa- 
tions {H — E)\'il)p^^) = with H = Hq + Vc and represented in the \k) space for the two 
protons. fTo is the S-wme Coulomb phase shift ctq = aigT{l+iri). The S'-wave amplitude 
As is given in terms of the effective range parameters as 



"n^N --^ + \roP^ H amNh{rj) -ipC^' 



(17) 



Figure 3: Dressed di-baryon propagator without Coulomb interaction (double line with 
a filled circle) at leading order. A single fine stands for the nucleon, while a double line 
represents the bare di-baryon. 




Figure 4: Diagram for the S'-wave NN amplitude without Coulomb interaction at leading 
order. The double line with a filled circle represents the dressed di-baryon propagator 
obtained in Fig. 3. 



where ac is the scattering length, ro is the effective range, and the ellipsis represents the 
higher order effective range corrections. Now it is easy to match the parameters ag and 
Us with the effective range parameters. Thus we have cr^ = — 1 and 

= ±J_f^ D ( ) - I (18) 

mjv V ^0 ' ' ^ ^ ~ I'^oP^ + amNh{r]) +ipC^' 

In Fig. 3, LO diagrams for the final deuteron channel, i.e., the dressed ^5*1 channel 
di-baryon propagators are depicted. Since insertion of a two-nucleon one-loop diagram 
does not alter the order of the diagram, the two-nucleon bubbles should be summed up to 
the infinite order. Thus the inverse of the dressed di-baryon propagator for the deuteron 
channel in the CM frame reads 

iD;\p) = iat(E + At)+iy','^{zp) 



An 



An at At AnatE 

n- H o- + ZP 

rriNyt ruNyt 



(19) 



where we have used dimensional regularization for the loop integral and E is the total 
energy of the two nucleons, E ~ p^ jmf^. The dressed di-baryon propagators are renor- 
malized via the 5'-wave NN amplitudes. The amplitudes obtained from the diagram in 
Fig. 4 should satisfy 

4:7r % 

lAt = {-iyt){iDt{p)\{-iyt) ^ 4... „2 ' (2°) 
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(a) (b) (c) 

Figure 5: Diagrams for the pp fusion process, pp — > de'^Ve, up to NLO. 

where At is related to the 5'-wave NN scattering ^'-matrix via 

S-l^e^^s.-l^^^ = ^(^)A. (21) 
p cotOt — ip \ 2n J 

Here 5t is the phase shift for the ^5"! channel. Meanwhile, effective range expansion reads 

p cotSt = -7 + ^Pdil^ + p') + • • • . (22) 

Now, the above renormahzation condition allows us to relate the LECs to the effective- 
range expansion parameters. For the deuteron channel, one has at — —1 and 

mN\ pd 2 7 + ^p - 2Pd(7 +P ) E + B 

where Zd is the wave function normalization factor of the deuteron at the pole E = —B, 
and the eUipsis in Eq. (23) denotes corrections that are finite or vanish at £■ = —B. Thus 
one has [18] 

Zd - . (24) 

This Zd is equal to the asymptotic S-state normalization constant. It is to be noted that 
the order of the LECs yt is now of Q^^'^, and the deuteron state is also described by the 
renormalized dressed di-baryon propagator. 

4. Amplitude up to NLO 

Diagrams for the pp fusion process up to NLO are shown in Fig. 5. In the limit p ^ 0, 
we have the amplitude from the diagrams in the figure as 

A = -^^d)-^mGFVud9ATfi. (25) 

Here e^^^ is the spin polarization vector of the out-going deuteron, e(;) is the spatial part 
of the lepton current P in Eq. (1), and 



Tfi 



^7r7 CrfC 

- IPd 7^ 



- acJxHx) + kac{ro + pdW + 

4 2gAmp 



(26) 
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where 

Hx) = - - e^E^ix) , E,{x) ^ Tdt^, (27) 

with X = ctfTLp/'-f. We note that the amphtude Tjj vanishes at the p ^ hmit because of 
the overall factor C^. The approximation is taken by keeping p dependence in while 
ignoring higher order p/tun corrections in the remaining part. Since p/niN ~ 0.2 %, 
the contribution from the higher order pjm^ terms will be sub 1 % order, which can be 
neglected conservatively at the uncertainty level we are considering in the present work. 
Introducing a "standard reduced matrix element" [16], 



A(p) = 



'^/.(P)!, (28) 



we have a finite and analytic expression of the reduced matrix element A(p) in the p ^ 
limit as 

A(0) = ^ - acl [1 - xe^^E^] + 7«c(ro + Pd)l^ + ^^kA ■ (29) 

As mentioned above, we exactly reproduce the result of the effective range theory at LO, 
and have a higher order correction proportional to the LEG Iia at NLO in Eq. (29). 

5. Numerical results 

We obtain the matrix element A(0) in Eq. (29) in terms of the four effective range pa- 
rameters, ac-i To, 7 and pd, and the LEG Iia- The values of the effective range parameters 
are well known, but three of them are slightly different in the references. In this work, we 
take two sets of the values: one is ac = -7.8063 ± 0.0026 fm, ro = 2.794 ± 0.014 fm, and 
Pd = 1.760 ± 0.005 fm from Table VIII in Ref. [33]. The other is ac = -7.8149 ± 0.0029 
fm, ro = 2.769 ± 0.014 fm, and pa = 1.753 ± 0.008 fm from Table XIV in Ref. [34]. 
We take an average of numerical values of A(0) from the two sets of the parameters for 
our numerical result. The value of the LEG Iia should be fixed by experimental data, 
but there are no precise ones for the two-body system. We fix the value of the LEG Iia 
indirectly from the relative strength of the two-body matrix element to one-body one, 
52B = M2B/M1B = (0.86 ±0.05) % in Eq. (29) in Ref. [14]. This value has been obtained 
from the accurate potential model calculation for the two-body matrix element with the 
current operators derived from HBxPT up to N^LO where the two-body current oper- 
ator has been fixed from an accurate experimental datum, the tritium lifetime, for the 
three-body system. Thus we have 

= -0.50 ±0.03, (30) 

where we have used our LO amplitude as the one-body input. This is a good approxima- 
tion because the difference between the amplitude from the effective range theory, which 
is almost the same as our LO result, and that from accurate potential model calculations 
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Our result KR(NL0)[16] 


BC(N4LO)[17] 


Pot. model[ll] 


A\0) 


7.09±0.02 7.04~7.70 


6.71~7.03 


7.05~7.06 



Table 1: Estimated values of A^(0). The value in second column is our result. The values 
in third, fourth, and fifth column are estimated from the pionless EFT calculation up to 
NLO by Kong and Ravndal (KR) [16], that up to N^LO by Butler and Chen (EC) [17], 
and an accurate phenomenological potential model calculation [11], respectively. 

is tiny [12]. For other well known parameters, we use B — 2.224575 MeV, = 1.2695, 
rup = 938.272 MeV, and m„ = 939.565 MeV, and thus have 7 = 45.70 MeV, x = 0.1498, 
and Ei{x) = 1.465. 

Employing the values of the parameters mentioned above, we have Alo(O) = 2.641 
at LO, and Anloi(O) = 2.662 ± 0.002 from the first set of the parameter values and 
Anlo2(0) ~ 2.664 ±0.003 from the second one up to NLO. Thus we have an average value 

Anlo(O) = 2.663 ± 0.004 , (31) 

and A^lq(O) = 7.09 ± 0.02 where the estimated error bars mainly come from those of the 
effective ranges, tq and pa, and the LEG Iia- 

In Table 1, we compare our numerical result for A^(0) with those from other theoret- 
ical estimations, the pionless EFT without di-baryons up to NLO by Kong and Ravndal 
(KR) [16], that up to N^LO by Butler and Chen (BC) [17], and the accurate phenomeno- 
logical potential model calculation [11]. We find that our numerical result is in good 
agreement with the values from the former theoretical estimations within the accuracy 
less than 1 %. As discussed before, the uncertainties of the estimations from the pion- 
less EFT without di-baryon fields are still large, ~4.5 % for the KR's estimation up to 
NLO, and ~2.3 % for the BC's one up to N^LO, mainly because of the unfixed LEC Lia- 
Though the results in the previous pionless EFT calculations have the unfixed LEC Lia, 
we can directly compare our result of the amplitude A(0) in Eq. (29) to the expressions 
in Eq. (7) in Ref. [17], and fix the value of the LEC Lia- Assuming the higher order LEC 
KiA = 0, we have Lia = 1.27±0.12 fm^, which is consistent with our previous estimation, 
LiA = 1.18±0.11 fm^ in Ref. [21]. When comparing our result with that from the accurate 
phenomenological potential model calculation, we find that our result is overestimated by 
~0.5 % mainly because we have not included the important contribution from the vacuum 
polarization effect. 

As a last remark we would like to note that the precedent pionless EFT calculations 
include the higher order corrections in both wave functions and vertices with external 
probe. The contribution to A(0) from the wave functions read 2.51, 2.54 and 2.58 at 
LO, NLO and N^LO, respectively. In our calculation with di-baryon field, higher order 
corrections to the wave functions are incorporated naturally by the summation of effective 
range contribution to infinite order, which gives A(0) equal to 2.64. A great advantage 
of the pionless EFT with di-baryon field lies in that we don't need to care the higher 
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order contribution to the wave function, and it is sufficient to take into account only the 
corrections to the vertices with external probe. This advantage reduces the number of 
Feynman diagrams dramatically, and makes the calculation of higher order terms very 
simple. 

6. Discussion and conclusions 

In this work, we employed the pionless EFT with di-baryon fields including the Coulomb 
interaction, and calculated the analytic expression of the amplitude for the pp fusion 
process, pp —>■ de^Vei up to NLO. Employing the assumption to distinguish LO and 
NLO terms in the contact di-baryon-di-baryon-axial-current interaction, we reproduced 
the expression for the amplitude of the effective range theory at LO. The LEG Iia-, which 
appears in the contact di-baryon-di-baryon-axial-current interaction at NLO, is fixed by 
using the relative strength of the two-body amplitude to the one-body one, 823 1 which 
has been determined from the tritium lifetime in the HBxPT calculation, and thus we 
could perform the parameter-free-calculation for the pp fusion process. We find that our 
numerical result of squared reduced amplitude A^(0) is in good agreement with those of 
the recent theoretical calculations within the accuracy better than 1%. 

As mentioned in the Introduction, the current theoretical uncertainties for the pp fusion 
process is ~ 0.3% in the HBxPT calculaiton up to N^LO [14]. To improve our result to 
a few tenth % accuracy, it would be essential to include the higher order corrections in 
the modified counting rules discussed in the neutron beta decay calculation [30] : the next 
higher order corrections would be the a order and l/rriN corrections. It is known that the 
higher a order corrections, such as the vacuum polarization effect [9] and the radiative 
corrections from the one-body part [30]^, are significant, whereas the corrections from the 
1/miv terms would be Pc/rriN ~ 0.16%. It would be worth calculating the S factor for the 
pp fusion process in a few tenth % accuracy with the pionless EFT with di-baryon fields 
including those higher order corrections. 

Another issue that we would need to clarify is the value of the LEG Iia-, which has 
been fixed in this work by using the result from the HBxPT calculation. As discussed, 
e.g., in Refs. [14, 35], the LECs which appear in the two-di-baryon-axial-current or four- 
nucleon-axial-current contact interactions, denoted by Iia in the pionless EFT with di- 
baryon fields, LiA in the pionless EFT without di-baryon fields, and in HB^PT, are 
universal. In other words, those LEGs are shared by the processes, such as, the pp fusion 
process {pp de^Ug) [12, 13, 14, 16, 17], nn fusion process {nn de^Ue) [21], neutrino 
deuteron reactions {ued — > ppe~, v^d — > upv^) [36, 37], muon capture on the deuteron 
{^i~d — > nnUfj) [38, 39], radiative pion capture on the deuteron {n'd — > 71717 [40] and 
its crossed partner 76? — > nnn^ [41]), tritium beta decay [14], and hep process (p^He — > 
^Hee+z/e) [14]. If these LEGs arc determined by using the experimental data from one 
of the processes, the lattice simulation [42], or the renormalization group method [43], 
then we can predict the other processes in each of the formalisms without any unknown 

^Thc radiative corrections from the one-body part are quite significant, 2 ~ 3% level, and are con- 
ventionally included into the rcnormalized Fermi constant G'y ~ GpVud and the phase factor fpp in the 
estimation of the S factor for the pp fusion process. 
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parameters. In this respect, it may be worth fixing the LEG Iia in the same formahsm, 
the pionless EFT with di-baryon fields, from, e.g., the tritium hfetime extending our 
formahsm to the three-body systems with electroweak external probes. 
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